function [dm,pr_bachelors] = DO_combine_moments

    %----- ACS Educational Attainment -----%
    ACS = readtable('Raw Data Files/EduAttainment_2003.xls','VariableNamingRule','preserve');
    total = table2array(ACS(2,2));
    bachelors = table2array(ACS(2,9:12));
    frac_bachelors(1) = (sum(bachelors)/total);
    
    ACS = readtable('Raw Data Files/EduAttainment_2004.xls','VariableNamingRule','preserve');
    total = table2array(ACS(2,2));
    bachelors = table2array(ACS(2,9:12));
    frac_bachelors(2) = (sum(bachelors)/total);
    
    ACS = readtable('Raw Data Files/EduAttainment_2005.xls','VariableNamingRule','preserve');
    total = table2array(ACS(2,2));
    bachelors = table2array(ACS(2,9:12));
    frac_bachelors(3) = (sum(bachelors)/total);
    
    ACS = readtable('Raw Data Files/EduAttainment_2006.xls','VariableNamingRule','preserve');
    total = table2array(ACS(2,2));
    bachelors = table2array(ACS(2,9:12));
    frac_bachelors(4) = (sum(bachelors)/total);
    
    ACS = readtable('Raw Data Files/EduAttainment_2007.xls','VariableNamingRule','preserve');
    total = table2array(ACS(1,2));
    bachelors = table2array(ACS(1,8:11));
    frac_bachelors(5) = (sum(bachelors)/total);
    
    ACS = readtable('Raw Data Files/EduAttainment_2008.xls','VariableNamingRule','preserve');
    total = table2array(ACS(1,2));
    bachelors = table2array(ACS(1,8:11));
    frac_bachelors(6) = (sum(bachelors)/total);
    
    ACS = readtable('Raw Data Files/EduAttainment_2009.xls','VariableNamingRule','preserve');
    total = table2array(ACS(2,5));
    bachelors = table2array(ACS(2,11:14));
    frac_bachelors(7) = (sum(bachelors)/total);
    
    ACS = readtable('Raw Data Files/EduAttainment_2010.xls','VariableNamingRule','preserve');
    total = table2array(ACS(1,2));
    bachelors = table2array(ACS(1,8:11));
    frac_bachelors(8) = (sum(bachelors)/total);
    
    ACS = readtable('Raw Data Files/EduAttainment_2011.xls','VariableNamingRule','preserve');
    total = table2array(ACS(1,2));
    bachelors = table2array(ACS(1,8:11));
    frac_bachelors(9) = (sum(bachelors)/total);
    
    ACS = readtable('Raw Data Files/EduAttainment_2012.xls','VariableNamingRule','preserve');
    total = table2array(ACS(1,2));
    bachelors = table2array(ACS(1,8:11));
    frac_bachelors(10) = (sum(bachelors)/total);
    
    ACS = readtable('Raw Data Files/EduAttainment_2013.xlsx','VariableNamingRule','preserve');
    total = table2array(ACS(1,2));
    bachelors = table2array(ACS(1,8:11));
    frac_bachelors(11) = (sum(bachelors)/total);
    
    ACS = readtable('Raw Data Files/EduAttainment_2014.xlsx','VariableNamingRule','preserve');
    total = table2array(ACS(1,2));
    bachelors = table2array(ACS(1,8:11));
    frac_bachelors(12) = (sum(bachelors)/total);
    
    ACS = readtable('Raw Data Files/EduAttainment_2015.xlsx','VariableNamingRule','preserve');
    total = table2array(ACS(1,2));
    bachelors = table2array(ACS(1,8:11));
    frac_bachelors(13) = (sum(bachelors)/total);
    
    ACS = readtable('Raw Data Files/EduAttainment_2016.xlsx','VariableNamingRule','preserve');
    total = table2array(ACS(1,2));
    bachelors = table2array(ACS(1,8:11));
    frac_bachelors(14) = (sum(bachelors)/total);
    
    ACS = readtable('Raw Data Files/EduAttainment_2017.xlsx','VariableNamingRule','preserve');
    total = table2array(ACS(1,2));
    bachelors = table2array(ACS(1,8:11));
    frac_bachelors(15) = (sum(bachelors)/total);
    
    ACS = readtable('Raw Data Files/EduAttainment_2018.xlsx','VariableNamingRule','preserve');
    total = table2array(ACS(1,2));
    bachelors = table2array(ACS(1,8:11));
    frac_bachelors(16) = (sum(bachelors)/total);
    
    ACS = readtable('Raw Data Files/EduAttainment_2019.xlsx','VariableNamingRule','preserve');
    total = table2array(ACS(1,2));
    bachelors = table2array(ACS(1,8:11));
    frac_bachelors(17) = (sum(bachelors)/total);
    
    pr_bachelors = mean(frac_bachelors);
    
    %============================== Data Moments in Table 5 =====================%
    
    dat = readtable('Raw Data Files/m1_data.xlsx','Sheet','Data','VariableNamingRule','preserve');
    dat_2003_19 = table2array(dat(88:291,5));
    m1 = 1-((1-mean(dat_2003_19))^3);
    
    dat = readtable('Raw Data Files/m2_NROU.xlsx','Sheet','Quarterly','VariableNamingRule','preserve');
    m2 = mean(table2array(dat(:,2)));
    
    
    m6 = -4.400;  % From Krowlikowski (2017), as stated in the paragraph after Table 5 
                  % This is reported in row 4 of Table 2 in Krowlikowski (2017)
    m7 = 2.133;   % From Huckfeldt (2022), as stated in the paragraph after Table 5
                  % This is reported in row 12 of Table 8 of Huckfeldt (2022) 
    m4 = 1.366;   % Source: https://www.bls.gov/news.release/archives/conemp_06072018.htm (See Table 7)
    
    
    dat_nb = readtable('Raw Data Files/m9_LEU0252917300A.xlsx','Sheet','Annual','VariableNamingRule','preserve');
    dat_ba = readtable('Raw Data Files/m9_LEU0252918500A.xlsx','Sheet','Annual','VariableNamingRule','preserve');
    wage_nb = table2array(dat_nb(:,2));
    wage_ba = table2array(dat_ba(:,2));
    m9 = mean(wage_ba./wage_nb);


    dat_stata = readtable('DO_ASEC_results.xlsx');
    
    m8 = table2array(dat_stata(1,10))*100;
    m10 = table2array(dat_stata(1,11))*100;
    m3 = table2array(dat_stata(1,12))*100;
    
    monthly_EU = table2array(dat_stata(1,13));
    m5 = 1-(1-monthly_EU)^3;
    
    dm = [m1; m2; m3; m4; m5; m6; m7; m8; m9; m10];

end